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Abstract 

We analyze the correspondence of many-particle and mean-field dynamics for a Bose-Einstein 
condensate in an optical lattice. Representing many-particle quantum states by a classical phase 
space ensemble instead of one single mean-field trajectory and taking into account the quantiza- 
tion of the density by a modified integer Gross-Pitaevskii equation, it is possible to simulate the 
superfluid to Mott insulator transition and other phenomena purely classically. This approach can 
be easily extended to higher particle numbers and multidimensional lattices. Moreover it provides 
an excellent tool to classify true quantum features and to analyze the mean-field - many particle 
correspondence. 

PACS numbers: 03.75.Lm, 03.65.Sq 
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Introduction. Ultracold atoms in optical lattices provide an excellent model system for 
the study of fundamental condensed matter problems, since experimental parameters can 
be controlled in a wide range and with high accuracy. If the atomic interactions are weak, 
bosonic atoms form a Bose-Einstein condensate (BEC). In the superfluid phase all atoms 
occupy the same delocalized quantum state and can be described by a macroscopic wave 
function. The dynamics of the condensate wave function is given in a mean-field approach 
by the celebrated Gross-Pitaevskii equation (GPE). However, mean- field theory takes into 
account only expectation values and neglects quantum fluctuations. Hence it is widely be- 
lieved to fail when fluctuations are of importance - striking examples are number squeezing 
and quantum phase transitions In this Letter we propose a phase space ensemble 

approach based on the mean-field limit including quantum fluctuations. If the atomic inter- 
actions dominate the dynamics, the system undergoes a quantum phase transition from the 
superfluid (SF) to the Mott insulator (MI) phase even at zero temperature due to quantum 
fluctuations. Since its first experimental observation 21 this transition attracts more and 

n 

more attention (see [3j and references therein). The MI phase is charaterized by a finite gap 
in the excitation spectrum and a vanishing compressibility. Number fluctuations are frozen 
out and the long-range phase coherence is lost. In an overall confining potential this leads to 
the appearance of a distinct shell structure. This pinning of the on-site occupation number 



to integer values has only recently been observed experimentally 
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If the optical lattice is sufficiently deep, the dynamics of the atoms is described by the 
Bose-Hubbard model Hamiltonian [5] 

H = J2 j (- J + a)a j+1 ) + ejhj + jh] \ , (1) 

where hj is the bosonic annihilation operator at the j-th lattice site and hj = Ojdj. As the 
Bose-Hubbard model is a genuine many-particle problem, approximations are of exorbitant 
interest. Numerically exact calculations can only be done for very few particles since the 
Hilbert space of the many-particle quantum states grows exponentially both with the particle 
number and the size of the lattice. 

The mean-field approximation for a weakly interacting Bose gas is usually derived within 
a Bogoliubov approach, considering the expectation values (aj) = y^N^pj only. Starting from 
the Heisenberg equations of motion for the operators dj and neglecting quantum fluctuations 
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(ojOj) w (a^)(aj) 2 yields the discrete GPE 



iijjj = -J(ip j+ i + Vi-i) + UNtylty (2) 

for the classical field This ansatz is exact for vanishing interaction U = and can be 
compared to the Ehrenfest theorem in single particle quantum mechanics. 

In this Letter we will propose a phase space approach to the mean-field limit and a 
generalization of the GPE which enables us to describe squeezed states and the SF-MI- 
transition purely classically. In contrast to the established Bogoliubov approach the many- 
particle state is represented by an ensemble of mean-field trajectories, taking into account 
also the higher moments of the many particle quantum state. The quantization of the 
density fij which gives rise to the Mott shell structure is achieved by the introduction of an 
integer Gross-Pitaevskii equation (IGPE). The phase space IGPE dynamics can be easily 
implemented and efficiently simulated even for higher particle numbers and two- or three- 
dimensional lattices. Moreover we will show that it provides an excellent tool to classify true 
quantum features and analyze the mean-field - many-particle correspondence like the break 
down of the mean-field approximation at unstable classical fixed points 6]. 

Quantum phase space. The key idea of the phase space IGPE approach is to map the 
quantum many-particle state onto an ensemble of trajectories in classical phase space, whose 
dynamics is given by the IGPE. Instead of quantum expectation values we will consider time- 
dependent classical ensemble averages. The initial ensemble is distributed according to a 
quantum quasi-probability distribution to transfer the characteristics of the quantum state 
onto the classical phase space. The Husimi distribution provides a non-negative and bounded 
function on quantum phase space, which is hence particularly suited for a comparison to the 
classical phase space distribution. Any many-particle quantum state can be represented 
by the Husimi distribution, which is given by the projection onto coherent states |{cj}), 

Q w (c 1 ,... 1 c M ) = \({c j m\ 2 - (3) 

In general, the structure of quantum phase space is determined by the operator algebra 
generated by the commutation relations of dj and at. For the M-site Bose-Hubbard model, 
the relevant algebra is SU(M), which reflects the conservation of the number of particles. 
The generalized coherent states [?J 

1 f^M 



are parametrized by the amplitudes at the lattice sites Cj which are normalized as ^ ■ \cj\ 2 = 
1. These concepts are most easily understood for the special case of just two lattice sites, 
where one can introduce the operators L x = |(aja 2 + a^ai), L y = f (aJjdi — a|d 2 ) and 
), which form an angular momentum algebra SU(2) with quantum number 
£ = N/2 (sj. The Hamiltonian (CQ) then can be rewritten in the Bloch representation as H = 
—2JL X + UL 2 up to a constant term. Since the SU(2)-algebra is topological equivalent to the 
(Bloch-)sphere S 2 (8, <fi), the coherent states and the Husimi distribution can be parametrized 
by the polar angle 9 and the azimuth angle <fi via the identification c\ = ^/(l + cos#)/2 and 
C2 = e 1 ^^/ (1 — cos6 l )/2. The dynamics of the Husimi distribution Q is then given by [q| 

-4 m( ^ + C08( ^l) g ' (5) 

The integer GPE. The classical Hamiltonian function TC of the mean-field dynamics, 
idipj/dt = &H/dip*, is given by the expectation value of the many-particle Hamiltonian (CQ), 
again neglecting quantum fluctuations: 

H = £ . (-J(r j+ x^ + cc.) + e^f + ^l^l 4 ) • (6) 

However, the quantum fluctuations, in particular the variance of the density, are of fun- 
damental importance in the MI phase. In principle these fluctuations act like an effective 
macroscopic variable driving the system towards integer filling. This can be already under- 
stood in the case of two wells: One can show that in the limit of strong number squeezing the 
quantum term ~ d 2 Q/d9dcj) in Eq. ([5]) simplifies and the dynamics becomes exactly Liou- 
villian again [9], however with a modified Hamiltonian function including a term describing 
quantum pressure. Generalizing this exact result to an extended lattice we propose a modi- 
fication of the GPE including the additional energy A7i = U/2 ^2j{(n 2 ) — (nj) 2 ) due to the 
quantization of the number operator. This term is only relevant if the interaction energy 
dominates the dynamics. In this case the lattice wells decouple and the many-particle wave 
function is given by the Gutzwiller ansatz as a direct product 1^) = ® ■ ■ ■ <g> \^m), where 
the local wave functions are written in a Fock basis \^j) = Yl n ^ \ n j)- Assuming that the 
local density is given by the mean-field expectation value iVl^ l 2 , the minimal ansatz for the 
local quantum state reproducing this density is given by 

|^.) = y/i - dj \rrij) + ^/d~e il3 \rrij + 1) , (7) 
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FIG. 1: (Color online) Comparison of quantum and classical dynamics in phase space approaching 
the hyperbolic fixed point for UN = 10, J = 1 and N = 20 particles. Left: Evolution of the Husimi 
density. Right: Dynamics of an ensemble of 150 classical trajectories at times t = 0, 0.3, 0.6, 0.9. 



where rrij is the integer part of the local density, i.e. N\tpj 



m,- + dj with rrij e N. 



The variance of the density is then given by (n|) — (nj) 2 = dj — d?. This ansatz presumes 
the smallest possible variance and thus the minimal energy correction A7i with respect to 
the GPE Hamiltonian (Q. Including this correction, Hamilton's equations yield an integer 
Gross-Pitaevskii equation (IGPE) 



iipj = -J(i> j+ i + ipj-i) + ejipj + U(mj + l/2)ipj. 



(8) 



Bogoliubov vs. phase space approach. As a first illustrative example, we consider 
the dynamics of the two- mode system. In the mean-field limit, the quantum expectation 
value (L)/iV is approximated by the Bloch vector 
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l^l| 2 -|^2 



(9) 
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The dynamics of s is restricted to the (Bloch-) sphere, reflecting the conservation of normal- 
ization. In the non-interacting case UN = the system performs common Rabi oscillations 
between the two modes. The anti-symmetric fixed point bifurcates for UN > 2 J, form- 




FIG. 2: (Color online) Comparison of the many-particle expectation value (L(t))/N (blue line), 
the Bogoliubov approximation s(t) (red line) and the ensemble average (s(t)) c \ for 1000 classical 
trajectories (dashed green line). Parameters are chosen as in Fig. [TJ 

ing two novel elliptic and one hyperbolic fixed point. The two novel fixed points have a 
non- vanishing population difference s z ^ 0, leading to the 'self-trapping' effect 8 1 ]. 

At the hyperbolic fixed point the classical dynamics becomes unstable and the Bogoliubov 
approach breaks down 6|. This is mirrored in quantum phase space as shown in Fig. [TJ The 
initial coherent state at 9 = 0.8ir and = diffuses in the direction of the unstable manifold. 
Therefore the quantum expectation value (L(t))/N penetrates into the Bloch sphere in the 
vicinity of the hyperbolic fixed point (marked by + in Fig. [2]), whereas the single classical 
trajectory s(t) stays on the Bloch sphere. However, this is by no means a breakdown of 
the mean-field approximation in phase space. As demonstrated in Fig. [TJthe dynamics of a 
classical ensemble closely follows the dynamics of the Husimi distribution. The dynamics of 
the quantum expectation value (L)/iV is well reproduced by the ensemble average (s) c i as 
shown in Fig. [2j 

The phase space approach can describe certain features of the quantum dynamics, such 
as squeezing and the dynamics at the unstable classical fixed point. Other aspects of the 
dynamics are genuine quantum like interference and tunneling in phase space. On the 
opposite, the GPE exactly reproduces the dynamics of the expectation value (L) for U = 0. 
However, the system is not necessarily 'classical' (cf. [10)). 

The superfluid to Mott insulator transition. The SF-MI transition is considered 
to be not explicable within mean-field theory because it is driven by quantum fluctua- 
tions. However, the phase space IGPE method takes into account fluctuations as well as 
the quantization of the density hj. Therefore it provides a fully classical description of the 
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FIG. 3: (Color online) Dynamics of the SF-MI transition in a harmonic trap with u = 0.17, M = 40 
and N = 26 particles calculated from classical ensemble averages: Variation of the parameters (a), 
eigenvalues of the classical SPDM (b), density pkk (c-d) and its standard deviation cr(pfcfc)(e-f). 
Figures (d) and (f) show pkk and &{pkk) i n the Mott phase at time t = t c calculated with the IGPE 
(blue circles) and the usual GPE (red crosses). 



the SF-MI transition. We consider the dynamics of N = 26 particles in an optical lattice 
with a superimposed harmonic trap described by the Bose-Hubbard Hamiltonian (CQ) with 
€k = u 2 (k — k ) 2 /2. We assume that the lattice is initially loaded with a pure BEC in 
the Gross-Pitaevskii ground state ijjj- The depth of the lattice is then increased adiabati- 
cally, suppressing the tunneling between the lattice sites, and decreased back again. This is 
described by a time-dependent hopping matrix element J(t) = Jq[1 + y 4e _ ^ _ic - )2// ^]~ 1 with 
Jo = 20, A = 10 3 , t c = 200 and t s = 40, while the interaction strength U = 1 is kept con- 
stant. We calculate the time evolution of an ensemble of 200 mean-field trajectories, whose 
initial amplitudes c^{t = 0) are distributed according to the Husimi function (J3J). A classical 
approximation of the single-particle density matrix (SPDM) = (a>\.ai) is then provided 
by the ensemble average pki = (c* k ci) c \. Figure [3] shows some characteristic features of the 
MI-SF transition calculated with the phase space IGPE, reproducing the results obtained by 



full many-body calculations 
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12| . Indeed one observes the occurrence of Mott shells with 



integer filling, p^k € N, and small superfluid regions in between. The density fluctuations 
are strongly suppressed in the MI phase, whereas they are significantly stronger between the 
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FIG. 4: (Color online) Dynamics of the SF-MI transition in a harmonic trap: Magnitude of the 
classical SPDM \p M \ at t = 0,130,200,400 (a-d). 

Mott shells. In the SF phase one eigenvalue of the SPDM is close to the particle number, 
indicating that the many-particle state is a coherent state described by one single conden- 
sate wave function ipj. In contrast, the MI phase is characterized by many eigenvalues of 
the same magnitude. This state cannot be described by a single condensate wave function, 
but by a phase space distribution. Figure H] shows the magnitude of the SPDM itself. One 
clearly sees that the coherences are lost in the MI phase at t — t c and mostly restored at 
t = 2t c . The SF-MI transition is reversible to a large extent, which is demonstrated by 
increasing J(t) again for t > t c . 

Let us notice that the suppression of density fluctuations and the loss and revival of 
coherences is already introduced into mean-field theory by the phase space ensemble ap- 
proach. However, the usual GPE cannot reproduce the Mott shells and the energy gap 
of the particle-hole excitations, because the quantization of the density is neglected. This 
becomes most obvious in Fig. [3] (d) and (f). The GPE predicts a smooth Thomas- Fermi 
density profile and a uniform suppression of density fluctuations in the MI phase. 

The occurrence of a gap in the excitation spectrum in the MI phase can also be under- 
stood within the classical phase space approach. Low-energy phonons cannot be excited 
as such collective excitations are impossible if the spatial coherences are lost. The remain- 
ing excitations in the MI phase are density fluctuations, which are discrete in the IGPE 
description. 
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Conclusion and outlook. Summarizing, we map the many particle quantum state to 
a phase space ensemble obeying a modified GPE including the backreaction of the quantum 
fluctuations onto the order parameter. The example of the SF-MI phase transition proves the 
power of the phase space ensemble method, providing an enormous alleviation of numerical 
effort and an illustrative insight into the many-particle dynamics. For example, the depletion 
of a BEC at a classically unstable point p] can be fully understood using phase space 
ensembles. The atoms are not lost - they are just redistributed to other modes. The phase 
space interpretation gives rise to fundamental questions: Which properties of a BEC are 
essentially quantum? Further work will be devoted to extended three-dimensional lattices 
and to a deeper analysis of the classical limit of many-particle quantum dynamics. Which 
criteria determine whether a many-particle system can be simulated classically? Advanced 
algorithms for many-particle simulations exploit the weak entanglement in ID quantum 
systems [yl and not the classicality of the state. 

Support from the Studienstiftung des deutschen Volkes and the DFG (GRK 792) is grate- 
fully acknowledged. We thank M. Fleischhauer and J. R. Anglin for stimulating discussions. 
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